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Abstract 

We study Hamiltonian difference schemes for scalar conservation laws with monotone 
flux function and establish the existence of a three-parameter family of periodic travelling 
waves (wavetrains). The proof is based on an integral equation for the dual wave profile 
and employs constrained maximization as well as the invariance properties of a gradient 
flow. We also discuss the approximation of wavetrains and present some numerical results. 
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1 Introduction 

This paper is concerned with oscillatory patterns in the nonlinear lattice equation 

2 iij + - = o, iieK, t > 0, jez, (l) 

which is, up to an appropriate scaling, a centred difference scheme for the scalar conservation 
law 

d T u + d^'(u) = 0, t>0, (eK. (2) 
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Although the lattice cannot be used for the approximate computation of (non-smooth) so- 
lutions of ([2]), there exist several reasons why it seems worth investigating the dynamical 
properties of ([1]) in greater detail. First, both the lattice and the PDE have the same Hamil- 
tonian structure, and ([I]) can be regarded as a variational integrator for ([2]) with respect to the 
space discretization, see appendix[A} Studying ([1]) therefore allows to understand in which as- 
pects quasilinear Hamiltonian PDEs differ from their spatially discrete counterparts. Second, 
the lattice is, similar to Korteweg-de Vries-type (KdV) equations, a dispersive regularization 
of ([2]), and hence it generates dispersive shocks instead of Lax shocks. These dispersive 
shocks describe the fundamental mode of self-thermalization in dispersive Hamiltonian sys- 
tems but are well understood only for systems that are completely integrable. The ODE 
system (TjQ) provides a class of non-integrable examples which can be simulated effectively, for 
instance using variational integrators for the time discretization. Finally, discrete conserva- 
tion laws provide toy models for more complicate Hamiltonian lattices. Fermi-Pasta-Ulam 
(FPU) chains, for instance, are equivalent to difference schemes for the so called p-system, 
which is a nonlinear hyperbolic system of two conservation laws. 

In what follows we focus on a particular aspect of the lattice dynamics and investigate 
special coherent structures, namely wavetrains. These are periodic travelling wave solutions 
to ([1]) and therefore linked to nonlinear advance-delay-differential equations. 

We emphasize that all subsequent considerations require a centred difference operator in 
(JJJ as only this one gives rise to a Hamiltonian lattice. Other discrete scalar conservation 
laws, as for instance upwind schemes, have different dynamical properties and are excluded. 

1.1 Dispersive shocks and wavetrains 

The relation between (fT]) and ([2]) manifests under the hyperbolic scaling 

r = et, £ = ej, Uj(t) = u(et,ej), (3) 

where < e <C 1 is the scaling parameter, and r and £ denote the macroscopic time and par- 
ticle index, respectively. On a formal level we can expand differences in powers of differential 
operators, i.e., 

u(£ + e) - u(£ - e) = 2%u(£) + £e 2 0fu(O + 0(e 5 ), 

so that, to leading order in e, the lattice dynamics is governed by the KdV- type PDE 

8 T u+ + ±e 2 df)$'(u) = 0, 

which is in fact a dispersive regularization of ([5]). 

As already mentioned, a key feature of all dispersive regularizations of ([2]) are disper- 
sive shocks. El El E] [25] p or illustration, and to motivate our analytical investigations, we now 
describe the formation of such dispersive shocks in numerical simulations of (fT]). For conve- 
nience we shall consider spatially periodic solutions with Uj = tij+jv for some iV S> 1, which 
then defines the natural scaling parameter e = 1/N. Starting with long- wave-length initial 
data 

«i(0) = u ini {j/N), 
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Figure 1: Numerical lattice solution with periodic boundary conditions for the data from 
(Q: Snapshots of uj against the macroscopic particle index j/N at several macroscopic times 
t; the vertical lines in the lower right picture mark the positions for the magnifications in 
Figure [2j Interpretation: Instead of Lax shocks the lattice generates dispersive shocks with 
strong microscopic oscillations. 



Zoom in at t=0.18, f=0.15 Zoom in at r=0.18, f=0.55 Zoom in at r=0.18, f=0.95 
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Figure 2: Magnifications of the oscillations from Figure [T] in three selected points: Black 
points are numerical data, Grey curves represent interpolating splines and are drawn for 
better illustration. Interpretation: The microscopic oscillations can be described by a modu- 
lated travelling wave. In other words, the local oscillations are generated by a single wavetrain 
whose parameters change on the macroscopic scale. 



where 'U'ini IS 3j smooth and 1-periodic function in £, we can solve ([1]) numerically by using 
the variational integrator (|28p . A typical example is depicted in Figure Q] for the data 

= \u 2 + \u\ iV = 4000, u ini (0 = l-isin(27rO. (4) 

For small macroscopic times r, we observe that the data remain in the long- wave-length 
regime, so we can expect that the lattice solutions converge as e — > in some strong sense to 
a smooth solution of • This convergence is not surprising in view of Strang's theoremESI j 
see also the discussion in Ref. |11}I17|. At time r ~ 0.04, however, the solution to the PDE 
forms a Lax shock, which then propagates with speed s according to the Rankine-Hugoniot 
jump condition s\[u]\ = |[<3?'(u)|. After the onset of the shock singularity, the lattice solutions 
do not converge anymore to a solution of ([2]), not even in a weak sense, but exhibit strong 
microscopic oscillations. These oscillations spread out in space and time, and constitute 
the dispersive shock. This phenomenon is fascinating from both the mathematical and the 
physical point of view because the formation of dispersive shocks can be regarded as a self- 
thermalization of the nonlinear Hamiltonian system. We refer to Ref. U [16] for more details 
including a thermodynamic discussion of dispersive FPU shocks. 
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The observation that certain difference schemes for hyperbolic conservation laws produce 
dispersive shocks is not new, see for instance Ref. |22|lll]. which investigate dispersive shocks 
in 

2ij + Vj(vj + i - Vj-i) = 0. 

This lattice also belongs to the class of equations considered here as it transforms via Uj = 
InVj into (pQ) with <£' = exp, which is the completely integrable Kac-von Moerbeke lattice^. 

The key observation in our context is that the oscillations within a dispersive shock exhibit 
the typical behaviour of a modulated oscillation. As illustrated in Figure [21 the local oscilla- 
tions in the vicinity of a given macroscopic point (r, £) resemble a periodic profile function, 
and we can expect this profile to be generated by a single wavetrain. The parameters of 
this wavetrain, however, depend on r and £, which can be seen from the fact that each 
magnification in Figure [2] displays a another amplitude, wave number, and average. 

The goal of this paper is to prove the existence of a three-parameter family of wavetrains 
for a huge class of nonlinear potentials <3?. Since the only crucial assumption we have to make 
is strict convexity of <£>, our results cover a large number of non-integrable variants of ([I]). 
Due to rigorous results for integrable systems, such as the KdV equation^ 1211 ano l the Toda 
chainEUE], we expect the parameter modulation within a dispersive shock to be governed by 
a variant of Whitham's modulation equations. 12QII2U The formal derivation and investigation 
of this Whitham system, however, is left for future research. We also do not justify that the 
oscillations within dispersive shocks take in fact the form of modulated wavetrains. For a 
rigorous proof in non-integrable systems we still lack the analytical tools; a reliable numerical 
justification would be possible, and was carried out for FPU in Ref. [I], but is beyond the 
scope of this paper. 

A travelling wave is a special solution to ([TJ with 

u J {t) = U{kj -ut), (5) 

where k is the wave number, ui is the frequency, and the profile U depends on the phase 
variable ip = kj — cot. In this paper we are solely interested in wavetrains, which have 
periodic U, but mention that one might also study solitons and fronts having homoclinic and 
heteroclinic profiles, respectively. 

Splitting U into its constant and zero average part via U(ip) = v + V(<f), we infer from 
([T]) that each wavetrain satisfies 

u-$fiV = V k *'(v + V), (6) 

where is defined by 

(V k P)(<p) = l[P(<p + k)-P(tp-k)). 
1.2 Main result and organisation of the paper 

Since the wavetrain equation ([6]) involves both advance and delay terms there is no notion 
of an initial value problem, and one has to use rather sophisticated methods to establish the 



4 



Snapshot at t=0 Snapshot at t=0.09 Snapshot at t=0. 14 




0. 0.5 1.0. 0.5 1. 0. 0.5 1. 



Figure 3: Numerical lattice solution for N = 4000, &(u) = u 2 and Umi{0 = sin(2-7r£). 
Interpretation: If changes sign, the lattices can generate modulated binary oscillation 
instead of dispersive shocks. 

existence of solutions. Possible candidates, which have proven to be powerful for other Hamil- 
tonian lattices, are rigorous perturbation arguments^, spatial dynamics with centre manifold 
reduction^ EH j critical point techniques^ ESIEZ]^ an( j constrained optimization^ ^1^3]. 
Our approach also exploits the underlying variational structure and restates Q as 

o-^'(Q) -arj = BQ, (7) 

where Q is the dual profile to be introduced in £ 12.21 Moreover, B is a compact and symmetric 
integral operator, and is the dual potential, i.e., the Legendre transform of This for- 
mulation allows to construct wavetrains as solutions to a constrained optimization problem, 
where a and —arj play the role of Lagrangian multipliers. 

The existence proof for wavetrains given below is based on a combination of varia- 
tional arguments (direct method) and dynamical concepts (invariant sets of flows). These 
ideas can also be applied to other Hamiltonian lattices and different types of coherent 
structures^ ESI 03] ^ Du t discrete scalar conservation laws are special since they are first order 
in time and real- valued. All other Hamiltonian lattices we are aware of are either second order 
in time or vector-valued, and allow for a simpler variational setting for travelling waves. 

Our main result guarantees the existence of a three-parameter family of wavetrains and 
can be summarized as follows. 

Main result. Suppose that $ is strictly convex and satisfies some regularity assumptions. 
Then, for each k £ (0, it) there exists a two-parameter family of 2ix -periodic wavetrains with 
DfR and frequency uj > such that the profile V is even and unimodal on [— ir, it]. 

The assumptions on <J? will be specified in Assumption [2] and the precise existence result 
is formulated in Theorem [8) Here we proceed with some comments concerning the choice of 
k and the convexity of <F 

1. Since ([5]) is invariant under (k, uj) (— k, — w), there is a similar result for k S (— ir, 0) 
with uj < 0. The cases k = and k = ±7r, however, are degenerate, see £12.11 

2. As time reversal in ([T]) corresponds to (k, uj) (k, —uj), the existence result can easily 
be adapted to the case of strictly concave <E>. 

3. Our proof does not cover potentials $ that change from convex to concave or vice versa, 
and it is not obvious what happens near zeros of Numerical simulations as presented 
in Figure [3l however, indicate that then a three-parameter families of wavetrains might 
not exist anymore. 
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The paper is organized as follows. In £|2]we summarize some elementary properties of wave- 
trains. We also derive the integral equation for the dual profile and discuss some normalization 
which allows to simplify the presentation. Sj3]is devoted to the proof of the existence result. 
To point out the key idea we first state an abstract result in Theorem [31 and show afterwards 
that the assertions we made are satisfied for wavetrains, see Theorem [8l Finally, in $3] we 
compute wavetrains by means of a discrete gradient flow on a constraint manifold. 

2 Preliminaries about discrete conservation laws 

In this section we summarize some elementary properties of wavetrains, derive the integral 
equation for the dual profile, and introduce some normalization. 

2.1 Elementary properties and special solutions 

First we observe that the wavetrain equation ([6]) is invariant under shifts p ~~» po, reflections 
<p ~~> —<p, and scalings 

k~^>k/s, u)~^*u/s, V(cp) ~* V(scp), v^v 

with s / 0. It is therefore sufficient to consider a fixed periodicity cell A and wave numbers 
k € A, so from now on we assume 

A = [—it, it]. 

Recall that a periodic function V : A — > R is even if V(<p) = V{—p) holds for all ip £ A, and 
that an even function on A is unimodal if it is monotone on [—it, 0]. 

Second, we notice that ([U]) degenerates for both k = and k = ±tt. In fact, wavetrains 
for k = are globally constant with Uj(t) = ^i(O), while wavetrains for k = tt are stationary 
binary oscillations with U2j(t) = uo(Q) and U2j+i(t) = Ui(0). Therefore, and since © is 
invariant under (k, u) (—k, —uj), we restrict all subsequent considerations to k € (0, it). 

For linear flux $' we can solve Q by Fourier transform. Specifically, linearizing ([6]) around 
v and restricting to unimodal and even profiles we find that uj depends on k and v via the 
dispersion relation 

co = v) = $"(v) sinfc, 

and that the unique profile is V((p) = a cos ip, where the amplitude a G R is the third 
independent parameter. Harmonic wavetrains furthermore exemplify that simple profiles can 
generate rather complex patterns for the spatial oscillations in the lattice. This is illustrated 
in Figure [U and in Figure [9] for the nonlinear case. 

The second case in which we can solve ([6]) explicitly are wavetrains with wave number 
k = vr/2. 

Lemma 1. Each wavetrain for k = tt/2 satisfies the Hamiltonian ODE 

^V(ip) = +^ v>sym {V(cp)), u^V(<p) = -& v>sym {V(ip)), (8) 
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Figure 4: Snapshots of harmonic wavetrains for several values of k. Although the carrier 
profile is always U(ip) = cos</?, the plots Uj = U(kj) against j S Z exhibit rather different 
spatial patterns. 



where 

V(<p) = V(<p + it/2), $ V:Sym (x) = |($(« + x) + $(v - x)). 

In particular, if <3> is smooth and convex, then for each v there exists a unique one- 
parameter family of wavetrains with unimodal and even profiles that is parametrized by 
Q = *v tB ym(y(0)). 

Proof. By a direct computation we infer from (|UJ) that 

w-£jV((p) = \$'{v + V(<p)) ~ l$'[v + V{fp - vr)) = -u,-fcV(<p + vr). 
Combining this with J A V dip = we find 

V(ip + Tr) = -V(ip) VtpeR, 

and hence flSJ) . Finally, since flSJ) is a planar Hamiltonian ODE the remaining assertions follow 
from standard arguments. □ 

2.2 Integral equation for the dual profile and normalization 

We define the dual profile Q of a wavetrain by 

&(v + V) = q + Q, q =/$'(« + V) dip, (9) 
where § abbreviates the integral mean value, i.e., 

■fvd<p=&r v(<p) d^ 

J A J-TT 

The key ingredient to our variational existence proof is to restate (J6j) as an equation for Q. 
To this end we consider the Legendre transform ^ of <£, which is well-defined and strictly 
convex provided that 3> is strictly convex 1 --, and satisfies 

(C) = sup (Cv - , $' O = O $' = id. 
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Using ^ and Q we now find that Q is equivalent to 

w-fc*'(q + Q) = V k Q, 
and integration with respect to tp yields 

uN/'(q + Q)-v)=AhQ, v= i *'(g + Q) dtp. (10) 

J A 

Here, v appears as a constant of integration and the integral operator A k is defined by 

i rip+k 

(A k Q)(p) = - / Q(<p) dtp. (11) 

z Jip-k 

With (|10p we have derived the dual formulation of Q; the main mathematical difference 
between both formulations will be discussed at the end of £13.11 



For the existence proof in £j3] it is convenient to normalize (|10p in two steps. At first we 
incorporate the parameter q into the nonlinearity by considering the normalized dual potential 

*,(C) = *(g + C)-* / (?)C-*(g), (12) 

which satisfies 

*£(C) = *"(? + 0, = *,(o) = o. 

This transforms (11011 into 



- 77) = AQ, t/ = « - tf'(g) = / %(Q) dtp. 

J A 

The second normalization step is motivated by the harmonic case and the observation that 
the formula for the phase speed a depends on the value of k. In fact, the linearization of (jlpp 
around q gives the dispersion relation 

sin ( k 1 

Q(tp) = acostp, uj = Q(q, k) = ^ , v = ^'(q), 

and we conclude that a = uj/k for < k < tt/2 but a = uj/(tt — k) for 7r/2 < k < tt, 
see Figure [5j We also notice that Ak has slightly different properties for < k < 7r/2 and 
tt/2 < k < tt. In particular, by (jlip we have 

Ar-itQ = -AkTQ (13) 

for all Q with ^ A Q dtp = 0, where T denotes the shift operator 

TQ = Q(-+tt). 

To complete the normalization we now replace uj by a, and ^ by either Bk or B n - k , which 
are defined for < k < tt/2 by 
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Figure 5: Dispersion relation for the harmonic case and the different relations between 
wave number k, frequency u, and phase speed a. Thanks to the symmetry (k, ui, a) <~~~> 
(— k, —uj, a) it is sufficient to consider k € (0, u). 

Specifically, for 0</c<7r/2we restate the dual wavetrain equation (jfOp as 

a{%(Q) - V ) = B K Q, k = k, a = u/k, (14) 

whereas for tt/2 < k < ir we write 

a(%(Q)-r}) =B K Q, k = vr - k, a = ui/(it - k). (15) 

Notice that for k = tt/2 both formulations agree and encode the ODE solution from Lemma[TJ 
In particular, in this case we have the further symmetry Q = —TQ- 

3 Variational existence proof for wavetrains 

In this section we prove the existence of a three-parameter family of wavetrains by showing 
that the dual integral equations (|14[) and (|15p possess corresponding three-parameter families 
of solutions. To elucidate the key ideas we start with an abstract existence result in §3. 1 1 
and check the validity of its assumptions in £13.21 

From now on we rely on the following standing assumption. 

Assumption 2. ^ is twice continuously differentiable and strictly convex with 

< c < *"(C) < c < oo 

for all C € M and some constants c, c. 

We remark that Assumption ([2]) holds if and only if <£, the Legendre transform of ^, is 
twice continuously differentiable and strictly convex with 1/c < <£" < 1 /c. This follows from 
basic properties of the Legendre transform^. 

Throughout the remainder of this paper, L 2 denotes the usual Hilbert space of square- 
integrable and periodic functions on the periodicity cell A, where the dual pairing and the 
integral norm are normalized via 

(Qi, Q2) = / Qi Q-z dip, \\Q\\l = i Q 2 dip. 

J A J A 
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Moreover, we define 



n = |Q G L 2 : dv? = 0| . 

3.1 Variational approach and abstract existence result 

We now prove the existence of a one-parameter family of solutions (Q, a, rf) € Tl x M + x M 
to the abstract wavetrain equation ([7]). To this end we assume that the operator B : H —> H 
is compact and symmetric, and that \£ is normalized by 

(0) = tf'(0) = 0. (16) 

Our variational approach to ([7]) is based on the functionals 

HQ) = \i QBQdap= |(Q, BQ>, W(Q) = / *(Q) d^, (17) 

JA JA 
and the constrained optimization problem 

maximize T on Af 7 = {Q e H : W{Q) < 7}, (18) 

where 7 > is a free parameter. We readily justify that ([7]) is the corresponding Euler- 
Lagrange equation for a maximizer Q G L 2 , where a and fj, = —arj are Lagrangian multiplies 
for the constraints W(Q) < 7 and J A Q dip = 0, respectively. 

The existence of maximizers Q € % can be proven by the direct method using weak 
compactness arguments, see the first part in the proof of Theorem[3J However, in order to gain 
more qualitative information about shape of the profile function Q we refine the optimization 
problem (|18p . For this purpose we introduce the corresponding (negative) gradient flow, that 
is the 'H-valued ODE 

^-Q = BQ- a(Q)V(Q), V(Q) = tt'(Q) - 77(g). (19) 



Here r > is the flow time, and 



V(Q) = -f ^'(Q)dip, a(Q) 
J A 



(V(Q), B® 

|2 



\W(Q)\\ 



2 



are two dynamical multipliers which guarantee that £. Q dip = and W(Q) = const hold 
along each trajectory of (fT9|) . Notice that, by construction, each stationary point of (fT9|) 
solves ([7]) with a = o~(Q) and rj = rj(Q), and vice versa. For later use we also mention that 
the explicit Euler scheme to (|19p with time step r > is given by the mapping 

Q ^ 1 T (Q) = (1 - t)Q + tBQ - to-(Q)V(Q). (20) 



We are now able to formulate the refined existence result. 
Theorem 3. Let C C H be any positive cone that 
1. is convex and closed, 
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2. invariant under the gradient flow (|19j) . 

3. satisfies supg gC F{Q) > 0. 

Then, for each 7 > the functional T attains its maximum on C n 7V 7 and each maximizer 
Q solves ([7]) with multipliers r] = rj(Q) and a = o~(Q). Moreover, each maximizer Q satisfies 
a{Q) > and W{Q) = 7. 

In preparation for the proof of Theorem we draw the following conclusions from As- 
sumption [2j the normalization condition (|16p . and the postulated properties of B. 

Lemma 4. The following assertions are satisfied. 

1. J- : L 2 — > R is well defined and weakly continuous, 

2. W : L 2 — > R is u>e// defined, convex, continuous, and Gateaux- differ entiahle with deriva- 
tive dW{Q) = $'(<3). We a/so We 

IsllOlll < W(Q) < ic||Q|| 2 , Q) > (21) 

/or allO ^ Q e L 2 . 

5. For eac/i 7 > 0, the set 7V 7 C H is convex and weakly compact in L 2 . It is also star- 
shaped in the sense that for each 7^ Q £ % i/iere exists a unique A 7 (Q) > suc/i i/iai 
AQ E My for all < A < A 7 (Q). 

Proof. All assertions follow immediately from the definitions of J 7 , W, and A/" 7 , see ()17p and 

(Hed. □ 



Lemma 5. The gradient flow (|19|) is we// defined on 7i \ {0} and conserves W. Moreover, 
T increases strictly on each non- stationary trajectory. 

Proof. Suppose that Q is given with V(Q) = 0, that means *f>'(Q) = const. The strict 
convexity of \& implies Q = const and Q € H gives Q = 0. From this we infer that o~(Q) is 
well defined for Q ^ 0. Consequently, and since the right hand side in (|19p is locally Lipschitz 
in Q, the ODE (|19|) is well posed in % \ {0}. By straight forward computations we now verify 
£W(Q) = and 

£T(Q) = (BQ - o-(Q)V(Q), BQ) = \\BQ\&- ( ^J^f > 0. 

In particular, we have -£pJ-(Q) = if and only if BQ and V(Q) are collinear, i.e., if and only 
if Q is a stationary point of (|19|) . □ 

We now finish the proof of the abstract existence result. 

Proof of Theorem^ With respect to the weak topology in L 2 , the objective functional J- is 
continuous and the constraint set C n W 7 is compact. The existence of a maximizer Q thus 
follows from basic topological principles, and by assumption we have F(Q) > 0. By LemmaS] 
we find 



A 7 (Q)V(Q) =T{\ 1 {Q)Q) <T(Q), 
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which implies A 7 (Q) < 1. On the other hand, Q G A/" 7 yields A 7 (Q) > 1, and thus we find 
A 7 (Q1 = 1, which means Q G <9A/" 7 . Moreover, Lemma [5] ensures that Q is stationary under 
the gradient flow ([15]) . and hence a solution to (J7|). Finally, testing ([7]) with Q we find 

2JF(Q)=a(y'(Q), Q), 

and J-(Q) > combined with (|2ip implies <r > 0. □ 

To conclude this section we mention that there also exist variational characterizations of 
solutions V to ((HI), but these are less feasible than the integral equation for the dual profile. 
The analogue to Q is 

aB- 1 V = &(v + V)-0, 9 = j $'(v + V)d<p, 

which has a variational structure, but a direct treatment is difficult as B -1 is not continuous. 
One can get rid off B^ 1 by making the ansatz V = BS with S G H and working with the 
functionals 

X{S) = \i SBS dcp, y(S) = i $(v + BS) dtp. 

J A J A 

However, the level sets of neither X nor y are weakly compact, and therefore it is not 
obvious how to set up a variational framework that allows to prove the existence of solutions. 
Finally, one might think about solving the equation V = B($'(v + V) — 9) /a by fixed point 
arguments, but then one has find a way to exclude the trivial solution V = 0. 



3.2 Proof of the main result 

Here we apply the abstract result from the previous section and prove the existence of a 
three-parameter family of wavetrains as claimed in the introduction. We therefore define 
C to be the positive cone of all L 2 -functions on A that are even, unimodal, and have zero 
average. This reads 

C = {q G % : Q(-f 2 ) = Q(f2) < Q(<pi) for almost all < < <f2 < tt}. (22) 
We next summarize some properties of the averaging operator Ak- 

Lemma 6. For all < k < tt the operator Ak maps T-L into itself and is symmetric and 
compact. Moreover, it maps C into C. 

Proof. The first three assertions follow directly from (|lip . It remains to prove that C is 
invariant under Ak for < k < tt/2; the claim for ir/2 < k < tt then follows from ()13p and 
the fact that T : C — > — C. Let Q G C be given, and notice that AkQ G H is even with 
weak derivative P = VkQ G L 2 . For <p G [—tt, —tt + k], the periodicity and the evenness of 
Q imply Q(ip — k) = Q(2tt + <p — k) = Q(—2tt — <p + k) where —tt < —2tt — <p + k < <p + k, 
so the monotonicity of Q in [—tt/2, 0] gives 2P(tp) = Q{ip + k) — Q(ip — k) = Q(ip + k) — 
Q(—2tt — <p + k) > 0. Similarly, for (p G [—tt + k, —k] and 92 G [—A;, 0] we have 2P(<p) = 
Q(ip + k) - Q(ip - k) > and 2P{tp) = Q(ip + k) - Q(ip - k) = Q{-ip - k) - Q{<p - k) > 0, 
respectively. In summary, AkQ is increasing on [—tt/2, 0], and since it is also even, the proof 
is complete. □ 
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Lemma [6] guarantees that both B K and B K are symmetric and compact, and map % — > H 
and C — > C. 

Lemma 7. For given geR and < k < n/2, Theorem^ holds with *$> q instead of ^ if C is 
given by ([22]) , and if B is replaced by either B K or B K . 

Proof. The cone C is convex and closed, and we have cos € C with (cos, B K cos) = 
(cos, £> K cos) = ^K _1 sin/t > 0. Therefore, it remains to show that C is invariant under the 
gradient flow (|19p . To this end we remark that C is invariant under the action of the operator 
Q i— > G(Q) — j'^G(Q) dip if and only if the function G : R — > R is non-decreasing. We now 
consider the Euler scheme of the gradient flow (|20p with time step r > 0, that is 

Q h-> X T (Q) = G TjCT (Q),7?(Q)(<3) + tB Q, 

where the family of nonlinear functions G Tt a , n : R — > R is given by 

G T , (7 , v (() = (-T(-Ta¥(0 + To-r 1 . 

Thanks to Assumption [21 for each bounded set U C R 2 there exists f > such that 
-d-G Tj a tV (() > holds true for all < r < r, all (eR, and all (a, n) G U. From this, 
the invariance of C under the action of tB, and the continuous dependence of a and rj on 
Q ^ 0, we conclude that for each ball B C L 2 \ {0} there exists r > such that X T maps 
B nC into C for all < r < f. The claimed invariance of the gradient flow now follows by 
passing to the limit r — > 0. □ 

We are now able to concretize the informal existence result from the introduction as 
follows. 

Theorem 8. Let < k < tt/2, q € R, and 7 > be given. Then, there exist profiles Q,Q E C 
along with multipliers a, a > and rj, rj such that 

a(%(Q)-rj)=B K Q, a(%(Q) - rj) = B K Q. 
In particular, the profile V = ty' q (Q) — rj^C solves 

with u = Ka > and v = rj + ^'(g), whereas V = ^' q {Q) — rj G C solves 

Z-feV = V*- K &(v + V) 
with u = Ka > and v = rj + W (q) . 

Proof. The existence of both Q and Q is a consequence of Lemma [7] and Theorem [3l The 
remaining assertions then follow by straight forward computations. □ 

We conclude this section with some remarks concerning the assumptions and assertions 
of Theorems [3] and [Sj 
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Figure 6: Numerical simulations for example (|23|) with several values of k. The plots show 
the graph of the dual profile, i.e., Q(y>) against <p G A. The lower right picture sketches the 
graphs of \I>" (Black), <&' (Dark Grey), and $ (Light Grey). 



1. The strict convexity of \I/ (or equivalently, of <D) is crucial for our existence proof. 
In fact, it is truly necessary for relating the integral equations (|14|) and (|15p to the 
wavetrain equation Q. It also plays an important role in the proof of Theorem [3] as 
it guarantees the weak compactness of 7V 7 . The assumptions about upper and lower 
bounds for ty" , however, are made for convenience and might be weakened for the price 
of more technical effort. 

2. It would be highly desirable to give uniqueness results that classify wavetrains up to 
phase shifts and scalings of the periodic cell. In analogy to the harmonic case, we 
conjecture that the three-parameter family from Theorem([8]) contains all wavetrains 
which satisfy the profile constraint V G C, but we are not able to prove this. Another 
open problem is the stability of wavetrains. 

3. Obviously, in Theorem [3] we can choose C = % to obtain solutions Q G A/" 7 to the 
original optimization problem (|18p . In this case, (|7|) is provided by the Lagrangian 
multiplier rule. If C, however, is a proper subset of H then the validity of © does 
not follow from the multiplier rule but is a consequence of the invariance properties of 
C. We also remark that numerical simulations indicate that each maximizer of J- in 
7V 7 is (up to phase shifts) unimodal and even, and hence contained in C. We therefore 
conjecture that (|18p and its refined variant from Theorem have the same solutions, 
but a rigorous proof is not available. 



4 Numerical simulations 

In this section we illustrate our analytical findings by numerical simulations, which are com- 
puted by the following implementation of the discrete gradient flow (|20p with *$> q instead of 
see (fT2|) . and B = Bk and B = Bk for k G [0, tt/2] and k G [vr/2, it], respectively: 

1. We sample the periodicity cell [—it, it] by grid points <p m = —ir + 2irm/M with m = 
1...M. 
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Figure 7: Numerical results for example (|24p . 
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Figure 8: Numerical results for example (|25p . 



2. We approximate Q by Q m = Q((p m ) and replace all integrals with respect to <p by 
Riemann sums. 

3. After each update step according to (f20]h we enforce the conservation of W(Q) by a 
scaling Q m (->■ XQ m , where the factor A is computed by a single Newton step. 

4. For given k and q, we initialize the scheme with Q m = acosip m , where the amplitude 
a determines 7 via 7 = ^ q (a cos ip) dip. 

Although this numerical scheme is rather simple it shows good convergence properties, pro- 
vided that M is sufficiently large and the time step r is sufficiently small. All simulation 
presented below are performed with r = 0.1 and M = 200. 

The first example concerns the data 

*"(0 = l + ^exp(-C 2 ), a = 5, q = 0, (23) 

where 7 is determined by a as described above. The numerical results are presented in 
Figure [6l which shows the dual profile Q as function of ip for several values of k. We also refer 
to Figure[9l which illustrates how the wavetrains appear in snapshots Uj = ^'(q + Q{kj — cot)) 
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Figure 9: Snapshots of uj against j = 1...50 for the simulations from Figure[6HH] a t a randomly 
chosen time t. Black points and grey curves represent the data on the lattice Z and the 
continuous carrier profiles U = v + V, respectively. Under the evolution all points move with 
constant speed along the grey curves. 



against j £ Z. A special feature of this example is that is even for q = 0, and this implies 
that ([7]) is invariant under Q ~» —TQ- We therefore find that the profiles are the same for k 
and 7T — k, and exhibit the further symmetry Q(ip + tt) = —Q(ip). 
The second example is computed with 

*"(C) = 1- Jarctan(2C), a = 3, q = 0, (24) 

and illustrates the generic case that ^ q is not even. In Figure [7] we therefore observe that 
the profiles do not satisfy Q = —TQ anymore, and are no longer the same for k and rr — k. 
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Finally, the third example uses 

vfr"(C) = iexp(2C)-iC-i « = 2, g = 0, (25) 

and is shown in Figure [7J Here is still convex but does not satisfy Assumption [2] due to 
^"(0) = 0. Nevertheless, the numerical scheme works very well and provides wavetrains for 
all values of k. 



A Lagrangian and Hamiltonian structures for scalar conser- 
vation laws 

For the readers convenience we summarize some basic fact about the Lagrangian and Hamil- 
tonian structures of discrete scalar conservation laws. Both structures are very similar to the 
respective structures for the KdV equation, compare for instance pQ, and originate from the 
skew-symmetry of centred difference operators. Specifically, for V with Vyj = ^(yj+i ~ Vj-l) 
the identity 

holds for all sequences y, y for which the series are well-defined. Notice that the one-sided 
difference operators V 1 * 1 with V^y^ = ±yj±i^fyj are not skew-symmetric but satisfy (V 1 * 1 )* = 

To derive the Lagrangian structure we assume that there exists y = yi{t) such that 
uj = Vy,-. This transform (TjQ) into 

Vy + W(Vy) = 0, (26) 
and reveals that the action integral is given by 



L(y, y) = -\f ^ Y, V^Vj dt ~ [ ^ E *( V %0 dt ' ( 27 ) 



where t^ n is a given final time. In fact, computing the variational derivatives d y L = — ^Vy 
and d y L = ^Vy + V3>'(Vy) we easily check that ([26]) is just the Euler-Lagrange equation 

_f td . L + dyL = . 

A special feature of L is that the canonical momenta d y L do not depend on y, and 
therefore we cannot apply the standard procedure to derive the Hamiltonian structure. There 
is, however, a non-canonical Hamiltonian structure. Using the Hamiltonian 

F(y) = J>(V yi ), 

iez 

the discrete scalar conservation law (pQ) can be written as 

Vy = d y H. 
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This is indeed a Hamiltonian equation, where V acts as non-canonical symplectic operator 
that corresponds to the symplectic product 



The PDE d^D has the same Lagrangian and Hamiltonian structures in the sense that all 
formulas from above remain valid provided that we (i) apply the scaling Q, (ii) replace 
sums over j by integrals with respect to £, and (Hi) employ the differential operator D = 
instead of V. In other words, the scalar conservation law ([2]) is the Euler-Lagrange equation 
to the action integral 

L(dy, y) = ~\ I d T yDy d£dr- / ° / <S>{Dy) d£dr. 

Jo Jm Jo Jr 

It is moreover equivalent to 

d T Dy + D$'(Dy) = 0, u= Dy, 

which has Hamiltonian H(y) = J R &(Dy)d{; and symplectic product (y, y') symp i = 
f^y'Dydt;. The analogies between the structures of lattice and PDE are - in view of the 
skew-symmetry of both V and D - not surprising and exemplify the more general principles 
laid out in [ID]. Moreover, if we replace D by the macroscopic difference operator 

we easily recover the scaled version of ([I]). In this sense the lattice ([!]) is just a variational 
integrator of the scalar conservation law ([2]), i.e., a semi-discrete scheme that respects the 
underlying variational and symplectic structures. 

The very same idea also allows to derive variational integrators for the time discretization. 
If we replace the continuous time derivative in (|27p by ^/i _1 (yfc+i, j — Vk—i,j)i where h is the 
time step size and the index k refers to the /c th time step, we obtain the discrete action 
integral 

L = -^^(vk+hj - Vk-i,j)(yk,j+i - ykj-i) -^2<s>{Uyk,j+i-ykj-i)), 

k,j k,j 

which depends on all the ykj's. Evaluating the Euler-Lagrange equation d Vkj L = 0, and 
replacing y by u, yields the fully discrete scheme 

Uk+i,j = u k -i,j - h$'(u k> j +1 ) + h&(u k ,j-x), (28) 
which also served to compute the numerical examples from £JU 
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